Refraction at Media with Negative Refractive Index 
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We show that an electromagnetic (EM) wave undergoes 
negative refraction at the interface between a positive and 
negative refractive index material. Finite difference time do- 
main (FDTD) simulations are used to study the time evo- 
lution of an EM wave as it hits the interface. The wave is 
trapped temporarily at the interface and after a long time, 
the wave front moves eventually in the negative direction. 
This explains why causality and speed of light are not vio- 
lated in spite of the negative refraction always present in a 
negative index material. 

PACS numbers: 78.20.Ci, 41.20.Jb,42.25.-p, 42.30.-d 

Veselago 1 predicted that lossless materials, which pos- 
sess simultaneously negative permittivity, e, and negative 
permeability, /i, would exhibit unusual properties such 
as negative index of refraction, n = —y/ep, antiparal- 
lel wave vector, k, and Poynting vector, S, antiparallel 
phase, v p , and group, v g , velocities, and time-averaged 
energy flux, (S) = (ii)v g , opposite to the time-averaged 
momentum density (p) = (w)k/tj, 2 where (u) is the time- 
averaged energy density. Furthermore if these materi- 
als are uniform, k, E, H form a left-handed set of vec- 
tors. Therefore, these materials are called left-handed 
materials (LHM) or negative index of refraction materi- 
als (NIM). The quantities, S, u, p, refer to the composite 
system consisting of EM field and material. As a result 
of k and S being antiparallel, the refraction of an EM 
wave at the interface between a positive n and a nega- 
tive n material would be at the "wrong" side relative to 
the normal (negative refraction). In addition, the optical 
length, J ndl, is negative in a LHM. 

Following Pendry's suggestions 3 ' 4 for specific struc- 
tures which can have both e e / / and /i e / / negative (over a 
range of frequencies), there have been numerous theoret- 
ical and experimental studies. 5-8 In particular, Markos 
and Soukoulis 9 have employed the transfer matrix tech- 
nique to calculate the transmission and reflection prop- 
erties of the structure suggested by Pendry 3 and realized 
experimentally by Smith et al. 6 ' 7 Subsequently, Smith 
et al. 10 proved that the data of Ref. [9] can be fitted by 
length independent and frequency dependent, e e //, and 
/i e //. They found that in a frequency region both e e ff 
and Heff were negative with negligible imaginary parts. 
In this negative region, n was found to be unambiguously 
negative. These unusual results [3-10] have raised objec- 
tions both to the interpretation of the experimental data 
and to the realizability of negative refraction. 11 ' 12 

In this letter we report numerical simulation results, 



which clarify some of the controversial issues, especially 
the negative refraction considered as violating causality 
and the speed of light. 12 Our numerical calculations were 
performed on a well understood realistic system, which 
is essentially inherently lossless, namely a 2D photonic 
crystal (PC). The dielectric constant e is modulated in 
space and both the permittivity e and permeability fj, are 
locally positive. Specifically, the photonic crystal consists 
of an hexagonal lattice of dielectric rods with dielectric 
constant, e = 12.96. The radius of the dielectric rods is r 
= 0.35a, where a is the lattice constant of the system. A 
frequency range exists for which the EM wave dispersion 
is "almost" isotropic. For that range the group veloc- 
ity, v g , is antiparallel to the crystal momentum k. 13 In 
particular for the frequency of our simulations the angle 
between the group velocity and the phase velocity pg 
varies between 176° and 180°. We have studied system- 
atically the cquifrcquency contours of the system 14 and 
have obtained the effective refractive index, n, as a func- 
tion of frequency. The absolute value of n, \n\ = c|k|/w, 
the phase velocity, v p = (c/|n|)k , where k = k/|k|, 
and the group velocity v g = daj/d\i are defined the usual 
way. Then, it can be easily shown that v g = v p /a, 
where a = 1 + (uj/n)(dn/du>) = 1 + c?ln |n|/dlna;; given 
this last relation is natural to identify the sign of n with 
the sign of a in order to associate negative n with an- 
tiparallel v p and v g . 15 Furthermore, it can be shown 
that (S) = (u)v g , where the symbol ( ) denotes averaged 
value over time and over the unit cell. Thus, v p and 
v g being antiparallel leads to almost all the peculiar and 
interesting properties associated with LHM. 

Essentially this means that for the PC system a fre- 
quency range exists for which the effective refractive in- 
dex is negative, frequency dispersion is almost isotropic 
and < S > k < 0, i.e. the PC behaves as a left-handed 
(LH) system. Consequently a wave hiting the PC inter- 
face for that frequency will undergo negative refraction 
for the same reason a wave undergoes negative refrac- 
tion when it hits the interface of a homogeneous medium 
with negative index n (the component of k along the 
normal to the interface reverses direction). In Fig.l we 
plot our results for the effective refractive index n for the 
PC system versus the dimensionless frequency /, where 
/ = Loa/2nc = a/X and A is the wavelength in air. Notice 
that the effective refractive index n passes continuously 
from negative to positive values as / increases. At this 
point we would like to alert the reader about other con- 
ditions under which negative refraction can occur in the 



1 



air PC interface. For example, for some particular k, 
light bends "the wrong way" at the PC interface as a 
result of the curvature of the equifrcquency surfaces 16 in 
spite of < S > k, and therefore the effective refrac- 
tive index-, being positive. Also in a PC system with 
< S > k > coupling to a higher order Bragg wave can 
lead to a negatively refracted beam. We would like to 
stress that in both the cases mentioned above, the PC 
is right-handed (RH). So the origin of the negative re- 
fraction in these cases is unrelated with that occurring 
on the interface between a RH and LH medium. Thus, 
these cases should not be confused with the case that we 
study in this paper where the PC is LH. 

In our simulations a finite extent line source was placed 
outside a slab of PC at an angle of 30° as shown in 
Fig. 2. The source starts gradually emitting at t = 
a monochromatic TE wave (the E field in the plane of 
incidence) of dimensionless frequency, / = 0.58, with a 
Gaussian amplitude parallel to the line source. Employ- 
ing a FDTD 17 ~ 19 technique, we follow the time and space 
evolution of the emitted EM waves as they reach the sur- 
face of the PC and they propagate eventually within the 
PC. The real space is discretizcd in a fine rectangular 
grid, (of a/31 and a/54 for the x and y axes, respec- 
tively), that stores the dielectric constant, and the elec- 
tric and magnetic field values. By use of a finite time 
step, (5i=0.0128a/c, the fields are recursively updated on 
every grid point. This algorithm numerically reproduces 
the propagation of the electromagnetic field in real space 
and time through the interface. In Fig. 2 we present 
the results for our FDTD simulation after 200 simula- 
tion steps. (The period, T = 2-k/uo corresponds to 135 
time steps, T ~ 1355*). 

We found, as expected for a negative index n, that 
the incident beam is eventually refracted in the negative 
direction (sin Opc < 0). However, the most interesting 
finding is that each ray does not refract in the final di- 
rection immediately upon hitting the surface of the PC. 
Instead, the whole wavefront is trapped in the surface re- 
gion for a relatively long time (of the order of a few tens of 
the wave period, T, in our simulations); and then, grad- 
ually after this transient time, the wave reorganizes itself 
and starts propagating in the negative direction as ex- 
pected from the steady state solution. Thus, the interface 
between the vacuum (positive n = 1) and the PC (nega- 
tive n) acts as a strong resonance scattering center which 
traps temporarily the wave before gradually re-emitting 
it. This time delay (which is much longer that the time 
difference, 2to ~ 400<5i, between the arrival at the inter- 
face of the inner and outer rays) explains satisfactorily 
the apparent paradox of the outer ray propagating much 
faster than the velocity of light. 12 

In Fig. 3, we present a time sequence of the amplitude 
of the magnetic field of the Gaussian beam undergoing 
reflection and refraction at the interface between vacuum 
(n = 1) and a negative index (n = —0.7) material. We 
would like to stress that despite the fact that the fronts 



are obscured by the Bloch modulation the refraction of 
the wave (after a transient time) is clearly in the negative 
direction. We have also calculated the Poynting vector 
(Fig. 4) which shows that (in the steady state) the en- 
ergy flows in the same negative direction as shown in Fig. 
3c. Since one can prove in general that (S) = (w)v g , it 
is clear that the group velocity in the PC is along the 
direction indicated in Fig. 3e. We present the results 
in terms of the time difference, 2to, between the arrival 
of the outer and the inner rays at the interface. Figure 
3a shows the results for t = 2M ^ 3.7T ~ 500<Si. No- 
tice that no refracted front has developed yet. Fig. 3b 
shows the results for t =4.5 to. Notice that the wave 
front has crossed the interface and seems to move along 
positive angles. Figure 3c shows the results for t =10.5 
to, where both the reflected and refracted wave fronts 
are shown. Notice there seem to appear two fronts for 
the refracted beam, one moving towards positive angles 
and one towards negative angles. Figure 3d clearly shows 
that the wavefront at t = 15io along positive angles has 
diminished and the wavefront along the negative angles is 
more pronounced. Also, notice the clear reflected wave- 
front. Finally, Fig. 3e for t =31 to clearly shows that the 
wavefront moves along negative angles, a behavior that 
continues for t > 31t - This is an unambiguous proof 
that negative refraction exists in media with a negative 
refraction index. This time sequence of the wavefronts 
shows that once the one end of the Gaussian wavefront 
hits the interface, it does not mean that the other side 
of the wavefront must move in "zero time" or "at infinite 
speed" to refract negatively (see Fig. lb of Ref. 12). 
There is a transient time of the order of 30t — 45T 
needed for the refracted wave to reorganize and move 
eventually along the negative directions. Notice that the 
wavefronts move finally along negative angles. The same 
is true for the directions perpendicular to the wavefronts, 
which can be considered parallel to the group velocity. 

We would like to note that Ziolkowski and Heyman 20 
have observed a similar transient time effect numerically 
in a ID model system with homogeneous negative per- 
mittivity e and permeability fi. In [20] it was shown that 
there exists a time-lag between the wave hitting the inter- 
face and the medium responding with negative effective 
index n indicating that causality is maintained. In the 
present paper we have shown that in a realistic struc- 
ture transient effect precedes the establishment of steady 
state propagation in the negative refracted direction. 21 

We note in passing that our time-dependent results 
are unrelated with the steady state solution obtained by 
Valanju et al. 12 In Ref. 12, a modulated plane wave, 
composed of two plane waves with different frequencies, 
is incident on the interface with a certain angle. Indeed, 
this steady state solution of propagation of a modulated 
plane wave, shows the interference fronts refracting in 
the positive direction, while the phase fronts refract in 
the negative direction. However Smith and Pendry 22 
have demonstrated that when two frequencies are super- 
imposed in the NIM with wavevectors that are not par- 
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allel - which is the case in [12]- propagation is crabwise 
meaning that the interference fronts are not the group 
fronts and should not be associated with the direction 
of propagation. In addition, Smith 23 has considered an 
incident modulated beam with a finite Gaussian profile. 
The steady state solution shows that the group velocity 
moves along negative angles, despite the distortion in the 
modulated fronts. Negative refraction is also obtained by 
Lu et. al. 24 where they consider more than two frequency 
components for the incident wave. 

In conclusion, we have shown by a FDTD simulation 
that an EM wave coming from a positive n region and 
hitting a plane interface of a negative n material re- 
fracts eventually in the negative angle direction. Our 
detailed time-dependent sequence shows that the wave is 
trapped initially at the interface, gradually reorganizes 
itself, and finally propagates along the negative angle di- 
rection. This transient time explains why negative re- 
fraction does not violate causality or light speed. 
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FIG. 1. The effective index of refraction, n, versus dimen- 
sionless frequency a/X for a 2d photonic crystal. 



3 




FIG. 2. An incident em wave is propagating along a 30' 
direction. The time is 200 simulation steps. 





FIG. 3. The magnetic field of the Gaussian beam undergo- 
ing reflection and refraction for (a) t = 2.5£o, (b) t = 4.5to, 
(c) t = 10.5to, (d) t — 15to, and (e) t = 31to. 2to is the time 
difference between the outer and the inner rays to reach the 
interface; to — 1.5T, where T is the period 2-k/lo. 
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FIG. 4. The magnitude of the Poynting vector for an EM 
wave propagating along a 30° direction. The time is 6200 
simulation steps, and is the same as the one shown in Fig. 3e 
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